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Branching fractions and direct CP asymmetries of charmless B decay 
modes at the Tevatron 

M. Morello a 

On behalf of the CDF Collaboration 

a Scuola Normale Superiore and I. N.F.N di Pisa - Ed. C, Polo Fibonacci, Largo B. Pontecorvo, 3 - 56127 
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We present new CDF results on the branching fractions and time-integrated direct CP asymmetries for 5° and 
B° decay modes into pairs of charmless charged hadrons (pions or kaons). The data-set for this update amounts 
to 1 fb _1 of pp collisions at a center of mass energy 1.96 TeV. We report the first observation of the B® — > K~ir + 
mode and a measurement of its branching fraction and direct CP asymmetry. We also observe for the first time 
two charmless decays of the At-baryon: A° — > pir~ and A° — > pK~ . 



1. INTRODUCTION 

The decay modes of B mesons into pairs 
of charmless pseudo-scalar mesons are effective 
probes of the quark-mixing (CKM) matrix and 
sensitive to potential new physics effects. The 
large production cross section of B hadrons of all 
kinds at the Tevatron allows measuring such de- 
cays in new modes, which are important to sup- 
plement our understanding of B meson decays. 
The still unobserved B® K~tt + mode could 
be used to measure 7 pQ and its CP asymmetry 
could be a powerful model-independent test of the 
source of direct CP asymmetry in the B system 
2J. This may provide useful information to solve 
the current discrepancy between the asymmetries 
observed in the neutral Acp(B° — > K + tt~) and 
charged mode A CP (B+ -> K+ir°) [3]. 

The B" -> 77+ir- and B° -> K + K~ modes pro- 
ceed through annihilation and exchange topolo- 
gies, which are currently poorly known and a 
source of significant uncertainty in many the- 
oretical calculations [4|5j . A measurement of 
both modes would allow a determination of the 
strength of these diagrams [5] . 

CDF II is a multipurpose magnetic spectrome- 
ter surrounded by calorimeters and muon detec- 
tors [7J. A silicon micro-strip detector (SVXII) 
and a cylindrical drift chamber (COT) situated 



in a 1.4 T solenoidal magnetic field reconstruct 
charged particles in the pseudo-rapidity range 
|?7 < 1.0. The SVXII consists of five concen- 
tric layers of double-sided silicon detectors with 
radii between 2.5 and 10.6 cm, each providing 
a measurement with 15 /an resolution in the 
azimuthal (0) direction and 70 /im along the 
beam (z) direction. The COT has 96 measure- 
ment layers, between 40 and 137 cm in radius, 
organized into alternating axial and ±2° stereo 
"super-layers" . The transverse momentum res- 
olution is <t Pt /pt — 0.15%pr/(GeV/c) and the 
observed mass- widths are about 14 MeV/c 2 for 
J /if) — > n + fx~ decays, and about 9 MeV/c 2 for 
D° — > K~tt + decays. The specific energy loss 
by ionization (dE/dx) of charged particles in the 
COT can be measured from the amount of charge 
collected by each wire. 

Throughout this paper, C-conjugate modes are 
implied and branching fractions indicate CP- 
averages unless otherwise stated. 

2. DATA SAMPLE 

We analysed an integrated luminosity J Cdt ~ 
1 fb _1 sample of pairs of oppositely-charged par- 
ticles with p T > 2 GeV/c and +Pt(2) > 
5.5 GeV/c, used to form B9-* meson candidates. 
The trigger required also a transverse opening an- 
gle 20° < A(f> < 135° between the two tracks, 
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to reject background from particle pairs within 
the same jet and from back-to-back jets. In ad- 
dition, both charged particles were required to 
originate from a displaced vertex with a large im- 
pact parameter do (100 /im < dg < 1 mm), while 
the B® s j meson candidate was required to be pro- 
duced in the primary pp interaction (d (B) < 
140 /xm) and to have travelled a transverse dis- 
tance L xy (B) > 200 /im. 

In the offline analysis, an unbiased optimiza- 
tion procedure determined a tightened selection 
on track-pairs fit to a common decay vertex. 
We chose the selection cuts minimizing directly 
the expected uncertainty of the physics observ- 
ables to be measured (through several "pseudo- 
experiments"). We used just two different sets of 
cuts, respectively optimized to measure the CP 
asymmetry Acp(B° — > K + ir~) and the branch- 
ing fraction B(B° S — > K~n + ), since those two 
measurements are the main focus of the analysis. 
For the latter, the sensitivity for discovery and 
limit setting [8] was optimized rather than the 
statistical uncertainty on the particular parame- 
ter, since the mode had not yet been observed. It 
was verified that the former set of cuts is also ad- 
equate to measure other decay rates of the larger 
yield modes (B° -> tt+tt", B° ->■ K+K~), while 
the latter, tighter set of cuts, is well suited to 
measure the decay rates and CP asymmetries re- 
lated to rare modes (B° s — > tt + tt~ , B° — > K + K~ , 
A° — > pn~ and A£ — > pK~). 

In addition to tightening the trigger cuts, in 
the offline analysis the discriminating power of 
the B® s j meson isolation and of the information 
provided by the 3D reconstruction capability of 
the CDF tracking were used, allowing a great im- 
provement in the signal purity. Isolation is de- 
fined as 1(B) = p T (B)/[p T (B) + £ iP t(*)], in 
which the sum runs over every other track (not 
from the B meson) within a cone of unit radius 
in the rj — <f> space around the B9-, meson flight 
direction. By requiring 1(B) > 0.5, we reduced 
the background by a factor four while keeping 
almost 80% of signal. The 3D silicon tracking al- 
lowed multiple vertices to be resolved along the 
beam direction and the rejection of fake tracks, 
reducing the background by a factor of two, with 



only a small efficiency loss on signal. The re- 
sulting 7T7r-mass distributions (see Fig. [1]) show a 
clean signal of B9-, — > h + h ~ decays. In spite of a 
good mass resolution (ss 22MeV/c 2 ), the various 
B? a ) h + h ~ modes overlap into an unresolved 
mass peak. 
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h + h ~ candidates passing all selection require- 
ments optimized to measure B(B® —> K~ir + ), us- 
ing a pion mass assumption for both decay prod- 
ucts. Cumulative projections of the likelihood fit 
for each mode are overlaid. 



3. FIT TO THE B DECAY-MODE COM- 
POSITION 

The resolution in invariant mass and in par- 
ticle identification is not sufficient for separat- 
ing the individual decay modes on an event-by- 
event basis, therefore we performed an unbinned 
maximum likelihood fit, combining kinematic and 
particle identification information to statistically 
determine both the contribution of each mode, 
and the relative contributions to the CP asymme- 
tries. For the kinematic portion, we used three 
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loosely correlated observables to summarize the 
information carried by all possible values of in- 
variant mass of the B candidate, resulting from 
different mass assignments to the two outgoing 
particles [9]. They are: (a) the mass calcu- 
lated with the charged pion mass assignment to 
both particles; (b) the signed momentum imbal- 
ance a = (1 —pi/p2)qi, where p\ (p 2 ) is the lower 
(higher) of the particle momenta, and qi is the 
sign of the charge of the particle of momentum 
Pi; (c) the scalar sum of the particle momenta 
Ptot = Pi + Pi- Using these three variables, the 
mass of any particular mode M\ 2 can be written 
as: 
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where mj (m 2 ) is the mass of the lower (higher) 
momentum particle. For simplicity, Eq. |T]) is 
written as a function of pi and P2, but in the 
likelihood it was used as a function of a and ptot- 
The simulated average values of as a func- 
tion of a for the eight — * h + h ~ and —=> ph~ 
modes are shown in Fig. [2l 

Particle identification (PID) information is 
summarized by a single observable n\(2) for track 

l(2\ defined as dE / d *K2)-dE/dx 1{2) (it) . 
\ /' ucnneu as dB /dxi m (K)-dE/dx im W wneie 

dE / dxi(2){Tr) and dE / dx\< 2 \(K) are the expected 
dE/dxi{2) depositions for those particle assign- 
ments. With the chosen observables, the likeli- 
hood contribution of the i th event is written as: 
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where: 

Cf n = R(M wv \a,ptot)Pj(a,Ptot), 
C k A iD = A(M nw ;c 2 ,m )P A (a,p tot ), 
4 in = e ClM -P E (a,Ptot), 
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Figure 2. Average M„r versus a for simulated 
samples of B° (top) and A° (bottom) candidates. 
The corresponding plots for the B° are similar to 
B? but shifted for the mass difference. 
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(7) 



£ A(e) = E ^f (E) < (E) ^™(«i, K2 |a,Ptot).(8) 

The various terms of the likelihood functions are 
described below. 

The index j runs over the twelve distinguish- 
able B? s -, — > h + h ~ and A" — » modes, and 
/j are their fractions to be determined by the 
fit, together with the total background fraction 
b. The background is composed of two different 
kinds: combinatorial background and partially- 
reconstructed heavy flavor decays. The combina- 
torial background is composed of random pairs 
of charged particle, displaced from the beam- 
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line, accidentally satisfying the selection require- 
ments, while the latter, referred as "physics" 
background, is composed of multi-body 6-hadron 
decays (i.e. B9 , — > pir/pK) in which only two 
tracks are reconstructed. The indices A(E) label 
the physics (combinatorial) background quanti- 
ties. The fraction of the physics background is 
given by j a and it is a free parameter in the fit. 

Each likelihood term, both for signals and 
backgrounds, is factorized into three different 
contributions: a) the conditional probability dis- 
tribution of the invariant mass M W7r given a and 
ptot (for the background M M is assumed to be 
independent of momentum), b) the joint condi- 
tional probability of PID variables k±, k,i given 
at, ptot for a determined particles hypothesis, j in 
the case of signals (Fj) and /,m in the case of 
background (Fi >m ), and c) the joint probability 
distribution of momentum variables a and ptot 

(-Pj(A,E))- 

R(M„„\a,ptot) = R(M n „-Mj(a, ptot), at, ptot) 
is the mass resolution function of each mode j 
when the correct mass is assigned to both tracks. 
The average mass M.j(a,ptot) is the value of 
M„„ obtained from Eq. (fl]) by setting the ap- 
propriate particle masses for each decay mode 
j and, by making a simple variable change, we 
obtain R(M T „ - Mj(ot, Ptot), ot,p to t) = R(Mj - 
Mb°(b°,A°)> a ?Ptot), where Mj is the invariant 
mass computed with the correct mass assignment 
to both particles for each mode j. a and ptot ap- 
pear explicitly in the last equation since they are 
useful to parameterize the dependence of the mass 
resolution by the momenta. 

The mass distribution of the physics 
background is parameterized with an "Ar- 
gus function", defined by the notation 
A(M 7rjr ; C2,mo) [TO], convoluted with a Gaus- 
sian distribution centered at zero with a width, 
in this case, equal to the mass resolution, while 
the combinatorial background with an exponen- 
tial function. The background mass distribution 
was determined in the fit by varying the param- 
eters ci, C2 and mo in Eq. (|5l6p . The function 
Pj(A,E)(a!,ptot) was parameterized by a product 
of polynomial and exponential functions fitted to 
Monte Carlo samples produced by a detailed de- 
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Figure 3. Tagged D° — » K~n + decays from 
D* + — > D°tt + — » [K~tt + ]tt + : a verification of 
the mass line shape by performing a 1-D binned 
fit where the signal mass line shape is completely 
fixed from the model (see text). 



tector simulation for each mode j, instead for the 
background terms was obtained from the mass 
sidebands of data. 

The mass resolution function R was parame- 
terized using the detailed detector simulation. To 
take into account non-Gaussian tails due to the 
emission of photons in the final state, we included 
soft photon emission in the simulation, using re- 
cent QED calculations (TTJ. The quality of the 
mass resolution model was verified in the simula- 
tion using about 500k D° — > K~ir + decays (see 
Fig.[3|). The mass line-shape of the D° — > K ~ir + 
was fitted by fixing the signal shape from the 
model, and allowing to vary only the background 
function. Good agreement was obtained between 
data and simulation. In Eq. ((4|), the nominal B°, 
and A° masses measured by CDF [12] were 
used to reduce the systematic uncertainties re- 
lated to the knowledge of the global mass scale. 

A sample of 1.5M D*+ -> D°ir+ -> [K-n+}7r + 
decays, where the D° decay products are iden- 
tified by the charge of the D* + pion, was used 
to calibrate the dE/dx response over the track- 
ing volume and over time, and to determine the 
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Figure 4. Tagged D° — > if ir + decays from 
D*+ -> D°7r+ ->■ [if-7r+]7r+: distribution of 
dE/dx around the average pion response, for cal- 
ibration samples of kaons and pions. COT sam- 
ples the amount of ionization charge produced by 
a track by measuring the time-over-threshold (ns) 
of the pulse on each wire associated to the track. 



Fj{l,m)( K i, K 2 \a,ptot) functions in Eq. (|7|8|) . In a 
> 95% pure D° sample, we obtained \Ao sepa- 
ration between kaons and pions (see Fig. HJ. The 
PID background term in Eq. is similar to the 
signal terms, but allows for independent pion, 
kaon, proton, and electron components, which are 
free to vary independently for physics (combina- 
torial) background. In Eq. ([8]) the indices I and 
m run over the four possible particles e, 7r, K, p 
and the fractions of different kind of particles 
wf~^ ,u?m^ are free parameters in the fit. Muons 
are indistinguishable from pions with the avail- 
able dE/dx resolution. 

4. FIT RESULTS AND SYSTEMATICS 

We performed two separate fits: the first us- 
ing the cuts optimized to measure the direct 
Acp(B° — > K + tt~) and the second to measure 
£>(£?" — > K~tt + ). Significant signals are seen 
for B° -> 7T+7T-, B° -> K+n-, and 5° 
K + K~ , previously observed by CDF [13]. Three 



new rare modes were observed for the first time: 
B° s -> K-JT+, A° b pir- and A° -> pK~ while 
no evidence was obtained for B® — > 7r + 7r _ or 
B°^K+K-. 

To convert the yields returned from the fit into 
relative branching fractions, we applied correc- 
tions for efficiencies of trigger and offline selection 
requirements for different decay modes. The rel- 
ative efficiency corrections between modes do not 
exceed 20%. Most corrections were determined 
from the detailed detector simulation, with some 
exceptions which were measured using data. A 
momentum-averaged relative isolation efficiency 
between B° s and B° of 1.07±0.11 was determined 
from fully-reconstructed samples of B®^ J/tp4>i 
B°^ D-W+, B a ^ J/iPK* , and B Q ^ D-w+. 
The lower specific ionization of kaons with re- 
spect to pions in the drift chamber is responsi- 
ble for a ~ 5% lower efficiency to reconstruct a 
kaon. This effect was measured in a sample of 
D + — > K~tt + it + decays triggered on two tracks, 
using the unbiased third track. The only cor- 
rection needed for the direct CP asymmetries 
A C p(B° -> K+tt-) and A CP {B° -> R-TT+) was 
a < 0.6% shift due to the different probability 
for K + and K~ to interact with the tracker ma- 
terial. The measurement of this correction has 
been achieved using a sample of 1M of prompt 
D° — > K~n + decays reconstructed and selected 
using the same criteria as B?<. — > h + h~ de- 
cays. Assuming the Standard Model expectation 
of A CP (D° -> R-TT+) = 0, the difference be- 
tween the number of reconstructed D° — > K~ir + 
and D — > K + ir~ provides a measurement of the 
detector-induced asymmetry between K + it~ and 
K~tt + final states. Since the same fit technique 
developed for the B9, — > h + h~ decays was used, 
this measurement provides also a robust check on 
all possible charge asymmetry biases of the detec- 
tor and dE/dx parameterizations. 

The B° -> K+K- and B° s -> vr+Tr" modes 
required a special treatment, since they contain 
a superposition of the flavor eigenstates of the 
Bg. Their time evolution might differ from the 
one of the flavor-specific modes if the width dif- 
ference Ar s between the mass eigenstates is 
significant. The current result was derived under 
the assumption that both modes are dominated 
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Table 1 

Results on data sample optimized to measure Acp(B° — > K + ir~) (top) and B(B® — > A~7r + ) (bottom). 
Absolute branching fractions are normalized to the the world-average values B(B° — > A + 7r~) = (19.7 ± 
0.6) x 10~ 6 and / a = (10.4 ± 1.4)% and f d = (39.8 ± 1.0)% [3 . The first quoted uncertainty is statistical, 
the second is systematic. N s is the number of fitted events for each mode. For rare modes both systematic 
and statistical uncertainty on N s was quoted while for abundant modes only statistical one. For the A° 
B(A b ^pir ) wag measured. 



modes only the ratio 



13(10-°) 



Mode 



N 3 



Quantity 



Measurement 



B° -> K+tt- 

B° -> 7T+TT- 

B° -> if+K" 



4045 ± 84 
1121 ± 63 
1307 ± 64 



B(B"^K- 7r+)-B(-B°^K+7r~) 
.X-7T+) + S(S°^X+7r-) 



B(B°- 
S(B°->7r" t 



B(B0->K+-k~) 
f 3 B(B°->K+ K~) 



-0.086 ± 0.023 ± 0.009 
0.259 ± 0.017 ± 0.016 
0.324 ± 0.019 ± 0.041 



5.10 ± 0.33 ± 0.36 
24.4 ± 1.4 ± 4.6 



B" -> K- 



B° -» K+K- 



A° -» pK~ 



230 ± 34 ± 16 



26 ± 16 ± 14 

61 ± 25 ± 35 

156 ± 20 ± 11 
110 ± 18 ± 16 



f, B(B l > 
fd 



*K+w~) 



B{B °->K+ tt~ )-B(B° s -*K- tt + ) 
B(B'l~tK+TT-) + B(B^~*K-T! + ) 
f d r(B°^K- Tr+)-r(B°^K+TT- 

r("B"-»x+7r-)-r(so^/f-7r+ 

f s B(B°^n+n-) 
7d~ B{B u -*K+-n~) 

B(B°->K+ K~) 
B(B"->K+tv-) 

B(Al^pn-) 



0.066 ± 0.010 ± 0.010 
0.39 ± 0.15 ± 0.08 
-3.21 ± 1.60 ± 0.39 

0.007 ± 0.004 ± 0.005 

0.020 ± 0.008 ± 0.006 
0.66 ± 0.14 ± 0.08 



5.0 ± 0.75 ± 1.0 



0.53 ± 0.31 ± 0.40 
(< 1.36 @ 90% CL) 
0.39 ± 0.16 ± 0.12 
(< 0.7 © 90% CL) 



by the short-lived B® component, that T s = Td, 
and Ar s /r s = 0.12 ± 0.06 PJTS] . The latter 
uncertainty is included in estimating the overall 
systematic uncertainty. 

The dominant contributions to the systematic 
uncertainty are the statistical uncertainty on the 
isolation efficiency (B® modes), the uncertainty 
on the dE/dx calibration and parameterization, 
and the uncertainty of the combinatorial back- 
ground model. The first is the larger systematic 
of all measurements with a B® meson decay (ex- 
cept for Acp(B° — > K~ir + )). The second sys- 
tematic, due to dE/dx, is a large systematic of 
all measurements, although the parameterization 
of the dE/dx is very accurate. The fit of com- 
position is very sensitive to the PID information. 
The third one is due to the statistical uncertainty 
of the possible combinatorial background mod- 
els and it is a dominant systematic for the ob- 
servables of the rare modes. Smaller systematic 
uncertainties are assigned for trigger efficiencies, 
physics background shape, kinematics, B meson 



masses and lifetimes. 

5. RESULTS 

The relative branching fractions are listed in 
Table [TJ where fd and f s indicate the production 
fractions respectively of B° and B® from frag- 
mentation of a & quark in pp collisions. An upper 
limit is also quoted for modes in which no signifi- 
cant signal is observed [16]. We also list absolute 
results obtained by normalizing the data to the 
world-average of B(B° — > A" + 7r~) 3]. The con- 
tributions from the likelihood fit for each decay 
mode are shown in Fig. [TJ 

We report the first observation of three new 
rare charmless decays B® — > K~ir + , A° — > pir~ 
and A° — > pK~ with a significance respectively of 
8.2cr, 6. Oct and 11.5ct. The significance includes 
both statistical and systematic uncertainty. The 
statistical uncertainty to evaluate the significance 
was estimated using several pseudo-experiments 
with no contributions from rare signals. 

The branching fraction of the newly observed 
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mode B(B° -v K~Tr + ) = (5.0±0.75± 1.0) x 10~ 6 
is in agreement with the latest theoretical expec- 
tation [17] which is lower than the previous pre- 
dictions 4 18 . Wc measured for the first time in 
the B® meson system the direct CP asymmetry 
of A CP {B° -> K-7T+) = 0.39 ±0.15 ±0.08. This 
value favors a large CP violation in B® meson de- 
cays, conversely it is also compatible with zero. In 
Ref. [5] a robust test of the Standard Model or a 
probe of new physics is suggested by comparison 
of the direct CP asymmetries in B° s — > K~ir + 
and B° — » K + tt~ decays. Using HFAG in- 

put 3 we measure -f— '- — *=o : : = 

0.84 ± 0.42 ± 0.15 where V is the decay width, in 
agreement with the Standard Model expectation 
of unity. Assuming that the relationship above 
is unity and using as input the B(B® —>■ K~n + ) 
measured here, the world average for Aqp (B° — > 
K + ir~) and B(B° —> K + n~) [3], we can estimate 
the expected value for Acp(B® — > K~ir + ) w 0.37 
in agreement with our measurement. 

The branching fraction B(B° S -> K+K~) = 
(24.4 ± 1.4 ± 4.6) x 10~ 6 is in agreement with the 
latest theoretical expectation [19 20 and with the 
previous CDF measurement |13j . An improved 
systematic uncertainty is expected for the final 
analysis of the same sample. 

The results for the B° are in agreement with 
world average values [3J. The measurement 
A C p(B° -> K+n-) = -0.086 ± 0.023 ± 0.009 
is the world's second best measurement and the 
significance of the new world average Aq P c - (B° — > 
K+n-) = -0.095 ± 0.013 moved from 6a to 7a. 

The updated upper limits and the absolute 
branching fractions of the currently unobserved 
modes B° — » K + K~ and B® — ► n + ir~ have been 
reported. The rate B(B° -> K+K~) = (0.39 ± 
0.16±0.12) x 10 -6 has the same uncertainty of the 
current measurements [3 a while the B® — » 7r + 7r _ 
upper limit (already the world's best [13j ) was 
improved by a factor 1.3, approaching the expec- 
tations from recent calculations [5|21j . 

We also report the first observation of two new 
baryon charmless modes Ajj — > pir~ and A° — > 
pK~. Wc measured £(Ag — > p7r _ )/S(A° — * 
pAT~) = 0.66 ± 0.14 ± 0.08, in agreement with 
the expectations from [22]. 
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